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ae 2 NGRODUCREGN 


lel GENYRAL. - - Marine navigation is the process 
of directing the movements of a vessel from one point to 
mnovner. Historically, niloting, dead reckoning, and 
celestiul methods have been used for this purvose. The 
earliest form of marine navigation was piloting, which 
came into being when man puided his vescel with respect 
to familiar landmarks. Dead reckoning probably followed 
as he attempted to predict his future positions in venturing 
beyond sight of land. While it is generally accepted that 
celestial bodies were used to steer by from almost the 
beginning, celestial navigation, as it is known today, hd 
to wait until information regarding the motion of celestial 
bodies was A Gace a 

The advent of World War II brought about the develooment 
Of Vurious nevigation:1 aids and electronic ydvositioning 
systems which enabled the navigator to fix his position 
without regard to cloud cover or weather conditions. Since 
this time the trend in navigation esninment his been toward 
Gbeainine preater positional sccuracy, Convinous posi vie. 
fixing capabilities, and longer range onerations. These 
improvements are being achieved by the use of electromagnetic 
and acoustic signals, esonhisticated electronic and celestial 


methods, and inertial systema. Automstic comput:.tion has 





become uw necessity in the more soph seio meee systems which 
require long and tedious calculations. 

The navigation requirements of the Navy are aquite 
different from those of the ordinary user. an ideal system, 
other than a completely self contained inertial system, 
would provide world wide covernge, all weather overation, 
an unlimited number of users, availability of frequent 
fixes, immunity to interference, passive user overation, 
and high positional eccuracy. Certain military operations 
Meguire oositional accuracies of 0.1 nautieal mile or better. 

Present long range navigation systems and their 

[2] 
Smaracteristics are presented in Table 1.1, where lone 
range is defined as having a position fixing capability at 
A range of approximately 400 nautical miles or grester. Of 
tnese, the inertial and acoustic doppler systems are 
essentially sopnisticnated dead reckoning systems with 
potentisl world coveruge and position Indicavane Cap ap Twine. 
However, at the voresent time and state of the art, the use 
of a navigational satellite system appears as tne best 
solution to the requirements outlined above. 

1.2 CONCSPTION OF A DOPPILER SATSLLITE NAVIGATTON SYST Me 
The conceovt of a dopoler navigational satellite system was 


by-product of the observation of the Soviet Union's first 


artificial earth satellite, Sputnik I, in the autumn of 





TABLE 1.1 - LONG-RANGE MAVTGATION SVST RS 


System 
Celestial 


Star 
teacker 


Ceneol 


bcbg 


Necca 


Loran-A 


Loran-C 


Omega 


SuteLlite 


beertial 


HeOustic 
Dopoler 


4istimated 


Wseful Range 


World-wide 


World-wida 


700-1200nm 


400-600nm 


200-500nm 


700 -900nm 


1200-1500nm 


5000nm 


World-wide 


4-12 hours 


> OO eae G 
weter denth 


~stimatead BEOGUuonc ¥y 
mC CUracy Range 
onm - 
Osonm —s 
Cs = Oaare 190-194ke 
2n'-263kKC 
2-6° 200-415kc 
280 -S25\ke 
0O.25-1.0Onm 10M ele)'de 
0.5-3.0nm emetic 
0.2-0.5nm 90-110ke 
O.5-1 .Onm LO-L5ke 
OS yr nep eis! LO0-400me 
2.Onm/hr - 
velocity 


1% of distance 
traveled 





Syetem 


Semestial 


Swor Tracker 


Consol 


RDF 


DeCCA 


Loran-A 


Loran-C 


Omega 


aoe. Loe 


Inmervial 


Acoustic Doppler 


TABL™ 1.1L 


(Gent ) 


Onin sro a 





Inexpensive equinment, easy to use; 
limited to 2 fixes/d-y, aff3cted by 
clouds and bad westher which «leo 
es Tine cee cee hele Ocul | Dae ome OOo x. 


Same as above excent equinoment 
Very Cxpehieive. 


Requires only a narrow bdand receiver 
to establish LCP. 


es CoOSsicet ai 
Capable of coarse fix only. 


World-wide coverage. 
beecons. 


Wasy to use, insufficient coverage, 
limited to adoute 740 mile range ap 
night because of sky wave interference. 


Basy to use, insufficient coversge, 
subject to sky wave interference. 
Less accuracy with sky wave. 


Hasy to use, 
coverage. 
wave. 


accurnve, inenutficaens 
Less accuracy with sky 
mxpensive installation. 


World-wide coverage. Presently in 
Ndvencec’ development. Will use wtnree 
frequencies to yield both coarse and 
fine grids. 


World-wide coverage. Tess sophisticated 
systems for comercial use now under study. 


Self centuined, must be checked with 
outside reference. Accuracy decreases 
in a few hours sat rapid rate. Very 
exnensive, requires expert maintenance. 


Sophisticsted high accuracy systems 
uLilize pyro-compuss, Velleci ty fiaeuc, 
and automatic track-plotting equipment. 
Limited to shallow water. 





ion. dr. We Me sGuiervanc De. G. C. Weiffenbuch of the 
Applied Physics Laboratory of Johns Hopkins lWniversity 
Listened to the signal generated by Soutnik snd noted 
the doppler shift in the transmitted radio frecuency. The 
shift is given quontitutively by: 
Ga) Ge (e//00 
where Af is the doppler shift from the satellite transmission 
freauency, Pis the range rate, ec is the velocity of 
electromégnetic provagation, nd f, is the satellite 
tranemission guentenen Hence, the mexsuring of doppler 
shift is compzrable to mensuring the range rate. This is 
the basis for the d2velopment of vosition determination 
equations derived in Chapter 3. 

Guier and Weirfendach deduced, and lster demonstrated, 
that accurate meusurement of the doppler shirt pxttern 
would permit the determination of a sutellite orbit. After : 
the first analysis and use of the dopoler tracking technique, 
Pr, fF. T. McClure of the Applied Physics Laboravory, fe aie 


tee tie technique could be inverted na that the co ples 


curve received from a sutellite whose position was known 


by 


could be used to deternine the vosition of tne receiver 


=> 


on the surface of the earth. It seamed only lcrical to 
utilizs this deduction as 4 navigation system. 
L.3 THe OPTRATICONAL RUIN Gy. =) = ne a ee 


constantly seeking improvements in the fundamental accuracy 





PoeaivreollOnmamaa in the operational character of havige«.tion 
een. Realizing that tne dopoler tecnniuve orovides a 
powerful tcol for simole Out very accurate navirntion, the 
@evstooument of a highly accurate overntionsl navigation 
system using artificinl satellites wis made on> of the Havy's 
Meroe Space objectives. The Chief of Navil Overstiona 
established %n ovesrationnl reouirement for a navirstion:l 
sateLlite as followe: 
Nevelop #2 satellite system to provide 
accurate, all-westher, world-wide nexvigation 
ror nuval esurfsace ships, sircraft, and 
submarines. [4] 
This renuirement was panerate?d by tne need for 4 hishly 
accurate position determination required by modern weapons 
systems. The satellite system was not developed to replace 
any of the oresent systems of navigation. Instead, it wus 
develoved to provide another navigation sid for determining 
positions of extremely high accurucy anywhere on the surface 
of the earth. 

Accordingly, the Navy propesed to the Advancec Research 
Projects Agency (ARPA) that a project be undertaken to 
exoloit the doppler techniaue for navigational purooses. 

The program began in 1958 as a feasibility study under 


contract with the Applied Physics Laboratory. This program 


was designated “Project Transit”. 


©? 





eee ee OL es te OF ea OPERATIONAL SYSTEM 


ge AG. 
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NERAT; CONFIGURATION AND OPERATTON - - After 
M@feninitial concention it dic not tuke lone to determine 
the configuration of an operational dovoler satellite 
nivigation system. The basic system consists of the 
followings five grouvs of eyuioment: 

1. Satellites. 

ce Tracking stations. 

5» Comouting center. 

4. Injection station. 
De Navigational ecuipment. 

Tne system is designed to employ four sateilites in 
meeer or near polar orbits atv an altitude of 600 miles ena 
a period of 110 minutes. With orbital planes separated vy 
45 degrees s fix is available at Least every 110 minutes. 
If the satellite onasses within 10 degrees of the zenith, 
the interval is extended to 220 minutes since longitude 
Can not be uccurately determined in this situaticn. At 
higher latitudes tne frequency of fix availability increases 
Ceonmciderably. “yolilcal satellite orbits for an operaricuer 

5] 
system are illustrated in Figure 2.1. The satellite 
transmits two rigorously coherent freouencies controlled 


meom one same OcCillevor. this allows | £arsc order 


Correction to be made for Lonospheric Yrelreacticn. wm 








FIGURE 2.1 - TYPICAL OPERATIONAL SATELLITE ORBITS 





wddition the sutellite transmits orbital information an 
Pome signals syncenronized with WTl-2,. The time sirnals 
are obtuined from 9 clock which is regulated by the 
satellite oscillator. 


The overall plan of a typical system operation is 


At 
eves 


illustrated schematically in Figure 2.2. f navigation] 
Saveliite is shown at three different time periods during 
merit. During time ®eriod qT the satellite signxis are 
received by 4 ground tracking station as the satellite 
passes within radio line of sight. As each tracking station 


Ig 


40) 


ceives the transmitted signal, it determines the avopsrent 
doppler shift and reds out the satellite memory and time 
marker. These data are then transmitted in digital form by 
teletype to the comvuter center. 

After receiving a sufficient number of pass reports, 
the comouter center uses the dovvler daanta to compute an 
improved satellite orbit and then voredicts the satellite 
position for svery two minutes im the next 16 hour period. 
When this is comoleted, the computer generates the parameters 
of the orbit from which the position of the satellite at 
any point in time can be determined. The comouter also 
analyzes tne time errors to give the clock rate, and 
determines tne correction to both the clock setting and the 


clock rate. wlleof this wformavion isst0ien Sear cic 


Prey ecGlem Seal1 One 


(D 
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Nuring time period To, when the satellite is within 
range of the injection station, the satellite memory is 
erased, the new data is inserted into tne memory, and the 
clock is reset and regulated. The satellite immedi-tely 
retransmits this message to the injection station where it 
igs compared with the information originally transmitted. If 
the readback indicates that the injection was unsuccessful, 
the message is automatically retransmitted to the satellite. 
When accurate storage in the satellite memory is verified, 
Pace omemory is locked by a time clockwssAdditional informaticn 
is rejected until approximately twelve hours later when 
the satellite is again in optimum range of the injection 
station. 

After receiving the message from the injection station, 
the satellite continuously transmits the purameters of its 
orbit and its predicted orbital position at two minute 
intervals. During time period Ts the navigation equipment 
aboard a ship recosives the orbital varameters from the 
sovellite and neaconda the doppler shite for at Leocteratce 
2 minute time intervals. Comnuting equivment uses this 
information to calculate the ship's true latitude, lonpitude 


Guan aaner 


> 


2.2 ~FPRASTBRILITY DaReRMiINATION. = =euiae 123s eee. 


\ 


of the Transit Program was vrimarily concerned with proving 
the feasibility of such a dopnler satellite navigation 


system. This phase was under the direction of ARPA. 


Pa 





During this period four exverimental satellites ware 
designed. They were desipnated the 1A, 1R, 2A, and 2B 
imesvoecuvively. Seven pround tracking stations were nut 
into operation and others were in the process of being 
designed and constructed. An ete imenval computation 
center wag placed into operation and a time and frecuency 
standards laboratory was established at the Avplied Physics 
Laboratory. The requirements of the injection station 
were determined and the design of the station was undertaken. 
axtensive mathematical anslyses of computationsl procedures, 
lonospheric refraction, geodesy, orbit determinution, and 
position determination were initiated. 

On 17 September 19099 the Transit 1A satellite wag 
launched by means of a Thor-Able vehicle. The satellite 
in it's flight shroud is shown in Figure 2.3. Althougn 
the satellite failed to obtain orbit, sufficient data were 
obtained to verify the feasibility of satellite tracking 
und navigation by doppler analysis. However, adjectives 
concerned with the verification of the satellite design 
were not achieved. 

The Transit 15 satellite was similar to the 14 in 
appearance and function. The satellite consists of a 
shell divided into two hemispheres 36 inches in Giameter, 
a central supeort, amd “i Anstrtimenv tray wet oe scum. 


9 lamination of fiberglas with a honeycomb plastic filler. 
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iim rorvovides a strong, hignly insulated, nonmetallic 
meeucvure. ‘White paint on the shell surface hids in 
Demperatvure control. ‘The central support, also made 
Sue Laminated fiberplas, connects the hemispheres of the 
ehell and supvorts the aluminum instrument tray. The 
tray is tied to the shell by nylon lacing which absorbs 
the shocks produced during lnunching. A cutaway view 
of the Transit 1B satellite is given in Fipure 2.40 

The tray contains two transmitting systems, silver- 
Zine and nickel-cadmium batteries, two command receivers, 
and a de-svdin system. <A radiation and insulation shield 
is vlaced above and below the tray. Two banks of solar 
cells used for charging the nickel-cadmium batteries and 
two de-spin weights are mounted outside and around the 
Bouavor of theeshell. Broadband “Amtennas Tr a somal 
pattern are painted on the outside of the shell. 

The tray also contains an infrared scanner which 
was designed for snecial studies by the Naval Ordnance 
Test Center. A spin rate of approximately 180 rpm was 
jieerced Go the satellite imeem « teav. Alter 4 peace 
time the de-spin system stops the spinnine of the satellite, 
since spinning renderk doppler datu useless. 

The Transit 1B satellite was successfully launched 
on 13 April 1960 by means of a two stape Thor-Able-Star 
vehicle. The elliptical orbit hiswen inclinati cnet ons 


degr:es, 4 perigee of 255 miles, and un apogee of 470 miles. 
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Tie satellite oeriod is 95.5 minutes. s«lthougn the 

satellite was not injected into the olunned circular, 

2900 mile orbit, the ground tracking sand computation vroved 
4 


completely successful. Analysis of telemetered data 


meer caven that the satellite cenmonents functioned 


an) 


seeistactorily. All electrical and mechanical asvectea o 
tne satellite were proven including the solar cell oversutison, 
command off-on operation, de-snin devices, and temperature 
stability. Navigational experiments were conducted which 


guve position determinations to within + mile of first order 


aad a. 
Pest 
Peak ad 


geodetic control points. All calculations used in orbit 
and navigational fix determination employed tne first order 
Mmomospneric refraction correction obtained from comparing 
the two coherent transmitted frequencies. The accuracy 
obtained indicated the validity of this two freouency 
Cechniaque. Investigations of the earth's gravitational field 
and geodetic surface were also conducted. 

Whe Peansit TR savéeliauve ceneed radiating on Lia, 
1960. However, the results of this experiment achieved 
the goal of tne first phase of the Transit research program 
by conclusively oroving the system feasibility while 
demonstrating the ability to perform dopoler satellite 
navigation with an accuracy required for most military 
pom! CaAbLona. 

2.5 PROVOTYPS] DiV=@LOPmeNT. = = WRth Uilemerce esa 
of the Transit 1B satellite the Transit Program entered 


Nes 





the second onase of development in esrly May 1960 when 
the management of the Program wag formally transferred 
from ARPA to the Bureau of Naval Weavons. This phsse 
consisted of the development of an engineering prototype 
of the operational satellite nuvigation system. Consecuently, 
a series of research and develoomeant satellites were 
Constructed and pl»uced into orbit. The basic churacteristics 
and orbital information of these satellites is given in 
us 

Table 2.1. “ffort was continued in improvement of system 
design, simolification, and reliability. ‘Comnrehensive 
testing in environmental chamoders, and for vibration, shock, 
and acceleration were conducted on the satellite and it's 
components. -A comolete data storage and readcut system was 
developed including the design and conetruction of an 
injection station. Tracking procedurea and esuipment 
were improved élong with tne development of shipboard 
navigation equipment. Continued research was conducted in 
the aresxs of pravity, refraction, and geodesy resulting in 
improved accuracy. 

The Transit 2A satellite was placed into orbdit on 
ef June 1960. 1 Giffered Prom tne 1A and 1A in thai 
more soler calls, an imoroved telemetering system, and an 
electronic clock. The clock, runnine from the satellite 
oscillator, orovided the accurate time correlation Cerween 
Une brackine statirone “hich 1s reduired for acer - 
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TABLE 2.1 - TRANSIT i. WARGH AND DEVGTOPMENT 


as J ¥ 


So i Rey as LAUy Ne een, DIRING 1960 AND 1961 
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Nav 2A | Nav 3B 
JUNE 21, 1960 {| FEB 21, 1961°° 


Nav 4A | Nav 43 | TRAAC 
JUNE 29.1961 | NOV 15, 1961 | NOV 15, 1961°°° 





LAUNCH DATE: APRIL 13, 1960° 













“B” SYSTEM: 162-21éme 162-216me 162-216me 
FREQUENCIES “C° SYSTEM; 54-324me §4-324me 54-324me 54-324me §4-32émec 54-324me 
“1” SYSTEM: 150-400me 150-400me 






5 Ports In 10, u 
5 Porls In 10! 


1 Part in 10 : 
7 Parts In 10 









RIMS FREQUENCY NOISE AS “B” SYS. 162-216me | 5 Ports In 1019 
MEASURED AT TRACKING "C" SYS. $4-324me | 5 Parts In 10 
STATIONS: “2” SYS. 150-400me 
Delay Line 
Register 


MEMORY: TYPE 
: CAPACITY 3B4 Bits 2049 Bis 1344 Bus 


WEIGHT | 265ib, | 273.3, | 290.3 tb. | 175.1 ts 233i 
INITIAL PERIGEE ALTITUDE ee asim | 9sdim | 


INITIAL APOGEE ALTITUDE: 754 km 1070 km a Wk 


* Note: Nav 1B ceased radiating on July 11, 1960. 
** Nav 3B entered the atmosphere on March 30, 1961. 
*** Launched pickaback on Transit 4B. 





2 Ports In 1010 
2 Parts In 10 








Magnetic Shift Ferrite Cave 





None 
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je Praneit SR satellite, lLaunened on 21 Febrimgey i961, 
mcoroOraved on Additional electronic unit in the Germ of a 
commlete dats storage system civable of etoring a small 
numoer of bits of digital information. A storsge losding 
techniove and readout syratem was used to insure correct 
Gata insertion. This provided the first flight test of 
the data system including the injection station. 
Mnfortunately, tne 3R satellite only rems»ined in orbit s 
little over five weeks and the experimental results were 
quite limited. 

Tne Transit 44 and 4F satellites were launched on 


29 June 1961 and 15 November 1961 respectively. <A cutaway 
ner 
pee 


View or the Transit 4A satellite is given in Figure 2.5. 
The 4 series was extensively engineered to schieve an 
extremely stable oscillator, small size, snd increased 
Mmercaoility. Crystal oscillators of high uccuracy were 
carefully sped and selected. Their operationszl temperature 
was actively controlled by a smail henting coil contained 
in the thermal inertial element. This emnohasis on constant 
opersting temperature yielded oscillator stabilities 
ie bo One part in Toe 10” short invervals. Such #eceurac, 
gives tne satellite system a capability for extremely 
precise position determination. 

The internal ecuioment was miniaturized to reduce 


the satellite size and weight. This was to enable the 
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Mice Or tin inexpensive Scout launch vehicle. ‘The ouser 

miele Of solur cells is collé#psable and foids back “gainst 
Bie satellite outer shell assembly. This greavly reduces 
be eize of the sntellite ss shown in Figure 2.5. Once 

cwne satellite was put into oroit the outer shell of solar 
cells was to be mechanically extended. In practice this 

Was never achieved since the launch vehicle was always large 
enough to allow the satellite to be fully extended within 
me a shroud. The 4 series also contained an experimental 
radio isotope (RIPS) pvower source. 

In order to increase the satellite reliability, plated 
circuits and connection welding techniques were used through- 
out. In addition, entire sybsystems were imodedded in piastic 
tO introduce tne necessary ruggedness, strength, end 
eremocrabure inertia. 


All the experimental saatellites so far have had to 


o 


use isotropic antennuis because tney were not stabilized 
with respect to the earth. Since no antennn is truly 

isouropic the sipnal in somes directions is much wealge: 
than in otners. Ag the satellite rotates, it sometimes 


¥ 


sweeps 2 direrzction of low signnil strength across tne 
ground receiver and tne signsl lock can not be mrintained. 
This accounts for almost aif of the sulellite Dagcace. 
that can not be used. 

The Transit Resemrcn andssattitude Control aa) 


satellite wags luuncned pickuback oy tne same vehicle as 


Alt 





tne 4B satellite. Its primary objective was tecting a 
method of gravity attitude stxhilizgntion. It uses the 
er dient of the graivitutional field across the dimensions 
er the sgtellite in order te keep one fuce of the satellite 
paways towird the earth. This allows the use of a directional 
antenn: and Avoids the problem of low signal strength in 
the antenna pattern. 

As & rasult of the data gathered from the research 
eae development satellites the Dasic characteristics anc 
performance of a satellite navigation systom were validated 
and a final prototype operational satellite was develonoéecd. 
moe firat of the operational Vransit 5 geries was Launched 


Teaneatpeon Sater bare 


t 


mi l702, An arvlist™s cCon@enviom of th 
- 
ig illustrated in Figure 2.6. Details pertaining to the 
operational satellite and system are not readily avsilsbdle 
due to security restrictions. However, this satellite 
incorvorsates thse design features proven in tne previous 
experimental sstellites. Conspicious additions include 
solar cell pacdles and a anmping system. The solar cell 
paddles can be angularly oriented with respect to the sun 
and thereby insure a nearly constant output of vower. The 
damoing system is comvosed of 2 long boom with a small 
weight xttacned to a spring at the end of thea boom. The 


large displacements of the spring vrovides the necessary 


Gnumping of angular oscillations reouired to keep the 
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FIGURE 2.6 - PROTOTYPE OPSPATIONAL SATSLLTVS 
TYP aro 





Miueilitve gravity stabilized. 

Since July i962 the 7. &§. Naval Observatory has 
maintained equipment for the recovery of time signals from 
Transit type satellites to stuacy the problems asscciatesc 
fea providing continuous synchronization of "U-2.. Precise 
time is required for highly accurate orbit and position 
cetermination. 

Nh later model of the prototype operational satellite 
was designed to use nuclear energey as 4 main vower source. 
When this satellite was placed into orbit it was the first 
operational satellite to use this new energy in space. To 
date, a number of prototype satellites have been placed 
mao Oroit. | 

Bet THE OPSRATIONAL S¥Stal a= = nut y oi Voce 
three of the vorototyve satellites were incorporated into 
an operational system. A complete operational ground network 
was set up and is entir®2ly separate from the research and 
Gevelooment network. The researcn and Gevelopment network 
is still in existance anc is designed to provide bDasic 
information on tne size and ashave of the earth and it's 
gravity field. This information is periodically processed 
through a computer at the Avolied Physics Laboratory and 
occasionally new force terms are derived. Tne new terms 
are then used dy the overational syetm to imorove orbital 
updating and predictions. In July 1965 the name of the 


Transit System was officially chsng7d to the Navy Navigation 


at 





Satellite System. The use of Transit is no longer authorized 
within the Navy. <A discussion of the ovnerational system 

and it's comnonents is oresentecd in Chanter 4, following an 
investigation of the mathematical theory of satellite 


navigation presented in Chapter 5. 


wv 
i 





O« NAVIGATIONMESSAYTSLULIGTE THsOey 


Belt Gaeeaal. =— = We will. sssume an ertvificial 
meee lite in a Known DOl4ar Of near polar OrDiiaeEansmi Ue ime 
on two extremely stable coherent frenuencies in tne 100 to 
400 megacycle range with added uniformly svaced time signals. 
At a receiving station an apparent change in the satellite 
cransmission freouency is SO tee Me is caused by the 
motion of the sitellite and the receiving station. The 
amount is pronvorti@onal to the,velocity of aporouch or 
recession. The exact amount of this shift depends on the 
exact location of the receiving station with respect to ta 

ss 

mun of the satellite. This Ps illustrated in Figures: am 
The navigation eauipment integrates tne heterodyne frequency, 
which is in the tens of kilocycles range, between the 
satellite transmission and a reference oscillator in the 
navigation equioment. “Ye mote from Figure 3.1 that the 


Fa) 
pt 


sign of the frequency shift che ngees from plus to minus at 

the center of passage. ' The beat note does not distinguish 
the sign of the freosuency difference. In order to get 
correct measurements it is reauired that the beat freauency 
does not pass through zero. Therefore, the local oscilixtor 
freauency is offset from the satellite oscillator by 80 varts 
in 10°. ‘This offset exceeds the maximum dopoler est £0 Ga 


that the beat nota never changes sign. WUocn raceipt of each 


time signal the number of cycles counted during the previous 
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OBSERVER "CLOSE TO SUBTRACK 
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OBSERVER AT INTERMEDIATS DISTANCE FROM SUBTRACK 


frequency 


+A 


ime 
-Af a 


OBSERVER DISTANT FROM SUBTRACK 


PIGURE 5.1 - HRPRPYCT OF OBSERVERS LONGITUDE DISPLACSMENT 
PROM SATELELI1Ys SUBTRACK. 





time interval is recorded. The navigation eauipment uses 
the higher of the two freouencies to d2termine count and the 
comonrison between ths two for na first order correction for 
Lonospheric refraction. It is then egssumed that the rrcorded 
Beycie count is refruction free. 

O.2 REPRACTION CCRRECTION. = = Accurate dofoler 


measurements have utilized the dispersive nature of the 


momastners to eliminate the refraction contribution which 


55 
is expressed as a series expansion in the inverse Yrrequency. 
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(Gee) me f(t.) += Os) = ae Et. : 
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c(t) F D(t) 


= <a ee 
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By receiving two coherent frenuencies transmitted from the 
satellite tne frenuencies may be substituted into (3.1L) and 
the first unknown coefficient 'A't determined at each time 
Vow CO Ey where © is the range rate and c the velocity of 
electromagnetic propagation. With the use of only two 
freauencies the resulting cycle count still contains all 
contributions that depend on tne higher powers of (1/2 
In operational systems two frequencies sare generally 
adequate for most precision culculations. the tracking 
station receiver uses analog methodsto automatically make 
the refraction correction before the dopnler data ure 
digitized. wm furtier discuscion of lonosunerie Tei wigs 


is oresented in Section 5.2. 
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ooo OSM enmines LON GP THe TNUeGkePsD CYCLs CHUNT. = = 


Tne Sieur € 


ellite Uransmits a time signal at 1, and at To. Tne 
tCranamissicn times from the sztellite to the r-ceaiver are 
1, 2nd to. If we denove: 
if, = satellite transmission freousncy. 
l — Pererence trequency. 


f. = obdserved satellite transmission freouency 
affected by Goopler snift, 


I¢ = integrated cycle count, 
we can derive the fcllowing relationships. 


The heterodyne frequency observed at the receiver is: 


<= ry. 


and the integrated cycle count is obtained from: 


ao “ ey 
(3.3) ie = Cae ce 


Somitting tne integral and evaluating the last term: 


aigtty 
— Geckos bao) ea eer | 


Ee 


Goes: ) 


tt 
C2 


P1+ hy 


Tne integral of (5.4) is notning more than the number of 
cycles transmitted by the satellite during the uniform 


cio ib erw a. Lo = Lye 


dO 





Therefore, we can say: 


my + 
ire & 
2 


meee =—- -¢ (f/m — Tp «a es es 7 98) | 
eo. 6 ) IC eer Py) Pp j'Tgr*a) (THe, ) 


Ry rearranging terms (3.6) can be written; 


e..7) 16] A fe fi As 
Wheres 

a =e 

. M9 - Ty 


Integrition allows the use of range differences (s3e 
souation 3.9) instead of range rate in compututions. 

Oe4 RELECION OF DISTANCES AND INTEGRATED CYCL4 0 Vite 
Woe basic distance formula states that distance 19 eonei ace 
Velocity times time. dence, the difference in transmission 


times from the sateliite is proportionul to the difference 


in range at T) and T,. Solving (3.7) for jay 
OE IN eS Ue 


(3.8) Lt = 
P. 


OL 
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CA 


$9) 


If we Let r enual tne ranyve to the satellite w 1 express 
eo. pane terma-ef distmnce. 


Cee te - IC) 
(5.9) a = 


With this basic eouation relating distance and integrat2d 
cycle count we can Gerive enustions for position 
Geterminution in a rectangular coordinate system. 
SO SCMATIONS FOR rOsITIONSOSTIRVINATION. = = 
If we denote: 
r= true distance to the satellite. 


R = computed distsnce to the sateilite from 
a DR position. 


Anat »4 = ssenmed shio coordinutes Imes riaae 
aes hand orthoronal coorlinute system. 
(obtained from some other navigation 
system) 
x?,y8,zS — : secellite coordinates in the same 
= Ge system. 
i = 0,1,2,-....-N - where W is enual to the number 
Of Lime SoOints. 
en 3 
fowl O} r =R-+dR 


Using assumed values for the variables as opoosed to the 
true Values in (o.9%) and @rtfecrentiatinge (6.9) and we owe, 


with respect to the variables we arrive at; 


Cc 
- : a Z = — (Aer Ae 
oe) a ee ee ~ | fice 
i & 
C Cc 
ers ae Paes 
. : (edie ee 
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Uxpreseing the R snd dR terms in terms of our orthorfonal 


coordinate «yetom: 


opal: Bee ee ex M4 (ve = y Sey) 
( ) ; ( i 7 "4S 77 XY 7d y ) 
Aa - Z - eye 
ue n = 
L " E : 
( 5 i) dRs gi a har eas X am OG ) oP, ( Ys ee Mg = ify \ oe 
val 


en i’ z)az_ 


i 


merere * x, \y, and Az account for navigstor motion bétween 


observations. Similiar enuations can be developed for bem 


—_ 


and GRs.y° By substitution of values obtained aodove into 


(3.11) Ar can be expresced as: 


i a oan , ‘ ; ar S U 4 A As | 
oe 4: ) P= ts = Rae e pals as x es iAx) AX, a 
x 
/ $05 2 i 7S en - r | 
Con Met isy)d¥, + (4; - 2) - inz)ea. | 
— 
ap — 
— (Koy = a x cA 
us 
ee 
S a ae oe fe as = = = ‘4 wi | 
rs Y (itl )Ay) ay -4 (254 Z ~ (iri)Az)az,, 


+ 
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Using this exoression for Ar in (%.11) and rearranging terms 
we can form an observation type egquation where v is the 
residual formed by the adjusted value of Ar minus the 


GoveVenam acetates leap isemre nm Lr. 
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imu Conele the cosificients of (5.14) as follows; 


a. = —eoetiicicent of dX. 

1 n 
Die COcihieient en av 
ae meet (eee Laon an 


mye equntion takes the form: 
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3.15 vo = a8.Gh 22 0D ae oe a7 Ae yer 
( ae n af ae al Lak cal rl 
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Hm this equation we assume our cycle count is errorless 
and therefore neglect the oIC term. Tne navigator may 
also sssume that the frequency of nis reference oscillator 
is unknown but constant during the short interval of 
observation. in tne eocuation this is treated as the unknown 
dAf. The sbvove eouation contains four unknowns (aX; qv 
AZ» df )- If we have four satellite observations we can 
solve for the unknowns and thus locate our position. 
Eauation (3.15) has been called an observation type 
equation instesd of an odservution equation oecsause the 
Pevacmal Is Lormed in bermipgot” 72 wiiecend of (1 erme in 
Section 3.7 tns actual observation eouation is derived 
wher? v is the residual formed by the adjusted value of IC 
minus the observed value of Ic. It will slso be shown that 
for navigationsl purposes Youation (5.15) may be Used" ace 


actual observation eouation in comoutntions. 





If we assume the earth to be aporoximated by an ellivsoid 


Seer evollution then: 


fe 








(5.16) en ef 
ee 
ae ha 


We can reduce the number of unknowns to three by expressing 


aZ, in terms of dx, and qy: Dirrerentiating (5.16). 


\ 


i ba ee | 3] 
(o. 177) da, = > ee AA a ache | jie | 
7 
45 ae al i | 
The unknowns are now dX_, Cae, and dAf. Now only three 


observations of a satellite are reouired to fix a vosition. 
By substitution of the value of dZ, from (3.17) into (35.14), 
multiplying and gathering terms, we can form an cbhservation 
Mevecton eimilar to (5.15)-% 

il : 


— ice -~ X. - iax) + 


de dl 
Ri \ 


ono) Vv. = Revd e m= 





- (Y5..7¥,7 (itliéy)4 ee Zy- (it+l)f2) 
i+1 
ve | e Cc 
(l-e== a oy = — {7 hen © IC) - a cee 
Za ye f fe 
nj Ne: Yr 


This equation can be reduced to tne form: 


, = oe One 7} | ae 
(5.19) Vs a,dx 4 b.dY or oo L. 
where As eeerricilent of One 


Ds -- coefficient of on 


Co eeeciiicicnt ot Gia. 
C 
ao \m _ 
foe) a 1G, 
Be 6 


Green es 


imp praccice™ we constrain our position to bewon a sohere 


since we ure dealing with small quantities. 


Mouarion (5.0) "becomes: 





(3.20) ee ae 
nh n QO e 
and (3.17) becomes: 
Lee ni 
(5.21) a7. = a ae oe Cy 
n . n n 
a 2 
The maximum difference in the dZ, tarm resulting from 


this approximation cun be determined as follows: 


(ee) avAn = (1 23 jaz, 
(elLivsoid) (sphere) 
and 
2D 
AZ - aZ == = Geez 
i ig n 
(ellipsoid) (snher?) (sphere) 





Therefore, the difference in (3.17) due to constraint to 

4 ephere is -e7dZ "he term aZ. haa it! atost 
: . 19) 4 2 — ” n ( sphere) 6 ANU r7 mm B19) ‘ 99a ke gS erea QaS 
Peeue 1b the eougtor. Assuming that the navientor position 
ig known to an accuracy of one nautical mile (<: 1829 meters) 

e hd e se 

and using the satellite derived value of e’ from the Kaula 
mieripsolda of 1965, we find the maximum difference in dZ, 


BO 9A; 
AaZ, = (0.00669)(1829) or 212.2 meters. 


This result is considered negligible for the purpose of 
navigution. By using the sphericel approximation the 

(1 - oe”) terme in the coefficients a; and bs are eliminated 
Mead cde obsiervitLion equation retains the form of (oe 
This is the general astellite observation eountion. 

2-6 REDUCTION OF ORS@RVATION BOUNATIONS. = =) Tiree 
ooservations of 2 satellite will yield three observation 
equutions. In vractice these are used directly in the 
determination of the three unknowns. ‘“Vhen there #re more 
observations than unknowns the method of leust squares muy 
be apolied to find that system of vilues which is most 
probable from the observutions themselves. arbitrarily 
assigning a weight of unity to euch observation and adding 


the square sums of the coeificients we obtzein: 


(5.23) [vv] = [pal dXi+ 2 [2b] aXgdyy + 2 BE axyaAryte [AT] ax+ 
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[by] ayg+ 2 bel ay,aAr, + 2 [bt] oy, + fee] car? + 2 fer aar,+ [Lt] 


bes 


: ee ‘ls Doge 
Meeres (ud) = a,4,4+ 45,5955, ---°-+ 4,8 


By the method of leust sousres the most likely solution is 
when the sum of the squares orf the residuals VV, Inge 8 
mooimume. To obtain this condition we set the purtial 
derivatives of (3.23) with respect to the three variables 
equal to zero. This resvlts in the formation of three 


normal equations. 
[Aa] ax + [adj ay_+ [ac] df, + [at] = 0 
524) fab] dX, + fbb! ay_+ foe] cAt,+ jbL] = © 
fac] 6X + ibe] d¥_+ [ec] df, + fel! = 0 


@Zauations (3.24) are normal equations reduced from satellite 
ooservation equations. Regardless of the numbder of 
observation enuations we can slways reduce them to a numbder 
of normal equations equal to the numbder of unknowns. ‘The 
solution of the normal equations ee unique results for 
the unknowns dX,, AV» and dAL - Coie one values for dke 
and dY, we can solve for qa Adding theese corrections to 


tne origins] values for x? nd and Zy, Kives the rectangular 


ne 


coordinates of our adjusted position. 


(5.25) Ne eer ces 


14 
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The value obtained for dAf, may be used to to update the 
freouency of the reference oscillnrtor in the navigation 
equipment. 

It snould 6? again noted that the use of the least 
squares solution is udvoided in practice by using only 
three observations to solve directly for the unknowns. If 
more than three observations are obtuined the excess 
ooservitions may be neplected or tney may be combined to 
mmgoduce a total number of three. This is the most pracvical 
eolution since the advuntapes of the leist aquares solution 
are minimized by the small number of observations that ure 
aviilable during one satellite nass. During & méeximum pass 
time of approximately 15 minutes only seven observ:.tions 
can oe made. at the vrsesent time a decreuse in the time 
materval to gain more observations 1s now pricewenl @ue uve 
eomoutvutionzl time requirements and restricrioms ©1 tac 
gatellite, memory and transmission system. 

Sef DERI VAGBON OF HE OBSERVE rON ROUATICN. - => Gabe 
Section 3.5 an observation type equation was formed where 
v is the residunl formed by the sdjusted value of Ar minus 
the observed value of Ar. In the theoretical case the 


ObSerVvition ecilstren senouler oe cert von nie err is /Ol sume 


Be. 





observed auantity. Such an obsarvation eavation will now he 
develoovsed in terms of the intepruted cycle count where the 
residual is formed by the adjusted IC minus the observed IC. 
Tt will then be shown that for the purpvose of navigation 
the ectual observation equation can be replaced by “quation 
eel). 
Using the notation of Section 35.5 And denoting: 

IG. == adjusted integrated cycle count. 

IG, == observed integrated cycle connt. 


we can expresa IC as a function of the unknown parametera. 


(35.26) CO ee) 
and : 
fone) 10, = ret voor vr= Ic, - IG, 


The adjusted parameters can be expresses as 


== ie, 
X. Xi x 
WV ui) me 
ee ee 8 
A, > 4y \ 2 


AN ae : 
i, 7a 


where xX, y, 2, and F are corrections to the assumed navigator 
parameters. Assuming that the corrections are small comoared 
to the parameters we can exvand (5.27) using a Taylor's 
series whore terms of the second order ana sreater are 


neplected. 
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The form of the 2ctual observation equation for doppler 
satellite observations can be expressed as; 


Break er nae ae 
(3.29) \ Rs: X he yi -—=—2z+-—-F +10. = i¢ 


yx \Y \Z leet n D 
el weg! ea" Bl 





ieorder to evalwate the coefficients of the corrections we 


mous from (o.9)3 


and by substitution: 


i | 

r Q qQ 

Sek = | ae S _ Z ‘ a 2 

a 4 n 
S) a Zs = rS os 2 Sol oy 2 
jhe mes (xf - x)? (xP - ¥, 

i 
(2: = Zee a 


The coefficients ‘can now be determined from (3.30). 


\ ic Py 

= \ ae le : 

yn =e I a, Os 7 
n ‘ ie ow a 


1 : i 
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If we denote: 


at ae 

As =e OOetl LO memes “OL ax, 
% ante 

Ds = COST Tre Pont tO. GY 
PLA . # 

Ce = cocci ficient. om dZy 


we can write (3.29) in the form: 


| Sb aa, a Aes : 
Meet) 6 Vy = a OX 4b oY 4c d2t ATS At + 10) = ic, 
Cc 
If we miltiply (3.31) by a constant equal to ~ — we 
f 
arrive at: : 
C ce 
r Bee a : ee 5 ; a werd 
(3.32) a a, AX, + d,dY + 050%, , AraAr. 
r c 
Cc 
a (IC, = ie) 
r 
The coefficients of this equution sare equal to the 
Coefficients of (5.15) se om o.oo) Sivemeorein 
C 
Arn = (Rye, + 2) = — (APA, - 169) 
r 





ou ava 


C 
Ee 


(3.15) can now be expressed as: 


C 





C 
me?) Vy = 230K te Var CHAZ --— AreAr, - : 


P (IC, = IC.) 
Y r 


C 
Se Ae age ge) 
f is 
r 
From Nquatious (3.32) and (3.33) the relstion between the 
actual observstion equation and the observation type egqusetion 


based on terms of Ar can be expressed as: 





C C 
(3.354) - vo SS vo 4 + Geel yeaa) 
a 5h at n b 
“> ik 
\rele 
5259) ee Kv KAT (Af =e Ge) 
Dele | =e i n b 
Cc 
where K is a constant equal to -— . 
ie 
Bouation (5.65) shows that Gee is related to vie by a 


constant except for the second term on the right side. It 
will now be shown thnt this second term can be neplected 


for the purpose of navisablom en comm ono) Wernece. 


Nt = i - we 
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ange: 


ee eae! “rn 
a i a re 
ae) eS ile: 
Yrom this it can be determined that the term (a - _—*) 
is equal to is - ee We will now determine the 


approximate maximum magnitude ofr the second term. Knowing 


that li has a drift rats of less than 2 varus in 10l0 per 


day and assuming that we use the solution of SAS a Prom 


orevious observations to correct f We Cat ee ei eat 


Tr? 


Poocoxcimavely 2 parle in 101°, Therefore, the term 


ol sion has an approximate maximum value of (2 24422¢)2 


ol, sing a value of 186,000 miles/sec 


Qe 2.8 Darts in 7 
for ce and 400 mc/see for f,, the approximate maximum value 


of the second term of (3.35) can be determined. 





ry ; 186,000 2.8 9 
OC OAS se: Ope ek Rie | (AL xO.) 
sa b 4x 108 109 


OY S2702006e miles ore 2) eer ere 


Since this value can be considered negligibdle for 


navigational purpoces we can write: 


BS IC 
oa } Ve == KvVy 


Therefore, if we set vic, IC to a minimum it follows that 
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T ts ey 
ae YS | also foes to a minimum. Hence we are able to 
Mowe. UOjein pace Of io051) in our computations. 


2.0 8A RECVANGULAR COORDIVAT OS. — = The adjuaved 
Seortiinate valwes in va rectangular systemeare of little 
Oracticsl use to the navigatior. Positional coordinates 
are desired in terms of latitude and longitude. While 
the calculations for coordinate determination and 
transformation are generally well known, they are included 
to orovide a comolete development of a position determination 
from donpler satellite observations. 

In calculating the positions of navipational fixes on 
the surface of the earth it is necessary to consider these 
voOints as lying vnon some mathematical surfsce. The actual 
shape of the earth's surface ia quite irrepyular and is not 
adapted to the purpose of comnutsation. For this reason it 
is neceseary to adopt an assumed figure which will not d>vart 
from the true surface by an amount sufficient to produce 
serious errors in our results. The figure generally adopted 
is an ellipsoid of revolution which is defined by a given 
flattening (f) and a piven semi-major axis (a). 

Whe position derivation given in the preceeding sections 
is with respect to an aroditrary ripgnt hand orthogonal 
coordinate system. We will now define thie system with 
resvect to the ellipsoid of revolution aszumed to revresent 
the figure of the earth. ‘The origin te au the Cenvcr ee ee 


ellipsoid which is also the assumec center of mass. Tne 
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Z axis coincides with the mean rotational axis and the X 
eeie Intersects the Grescnwich Meridisn at the equator. The 
Y ayis completes the right h.nd system. The coordinate 
system thus defined is illustrated in Fifure 3.3. 

The relutions between geodetic latitude snd longitude 
and the defined XYZ coordinate system, which can be derived 


from Figure 5.5, are given in the following formulas. 





2 COs Soa 
o= 
a cos (Q) sin \ é 
eo. 07 ) LY — oo 
W 
OL ae) sin@ 
Aa—_3.sere.,-”-—— 


q hd 
vd 


ver 


hore “| is cefiner 41s: 
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Boe A 
sin? @ )2 


2.9 DETPSRALRATIONSOr BATI TODS 2h GONGITUDS. -s- 
Using values for X,, Yoo and Z. from (3.25) and the relations 
given in (3.37) we are able to determine the latitude and 
longitude on the surface of the given reference ellipsoid. 


The value for longitude can be solved directly from (3.57). 


bier der | 
mee) A = fant —— 
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Mt Secison 


semi major arie 


semi minor axis B 2 





first eccentricity 
geodetic latitude 


geodetic longitude 


FIGURE 3.3 - wARTH RECTANGIITLAR GOORDTNATE SYST=M 
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Meee eOlLuLiLOm tor Latitude is a Little more involved. From 
Figure 5.4 we observe that: 


N sin 


4 


io. 39 ) tan@ = ine 32\5 
where N is the radius of curvature of the surface of the 
ellipsoid in a plane through the normal and at right angles 
to the meridian. This is also referred to as the radius 
@eecurvature in the prime vertical. A problem arises Chace 
feos a function of ) - (N is equal to a/V! ) To solve 
(5.59) we must eliminate @ from the right side of the 
equation. Going back to (3.37) we can express Z in the 


following manner: 


(3.40) Zo sine eNMo ein! 


and N sin@ = @ -+-N e@ sin® 
i 


Z+N ee” sin® Z N e© sin® 
(5.41) tan = 
(xe Y%)2 (X% Y¥%)2 Z 


Substituting N(1 - e“) sin@ for Z in the denominator 


of (53.41) we arrive at; 
(C5542) @ = Tan *( Sear ame.) 6 
(x2 y2)2 1 - e 


The values of latitude and longitude obtained from (3.38) 


and (3.42) are used by the navigator to fix his position 
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FIGURE 5.4 - MaRTDIAN SECTION IN BARTH COORDINATE 
OMS TSM 3 


49 





Ormmbne surface of the ellinsoid. 

Reversing this procedure the novigator crn eusily 
Barenlitve his position in esrth rectingular coordinetes 
ion values Of latitude and longitude by aonlying (oso7) 

OSel1QO INERTIAL COORDINATS SYST3M. - - The derivations 
to this point have assumed that the navigator is piven 
satellite coordinutes in the defined earth rectangular 
coordinate system. In practice the navigator must calculate 
these coordinates from given orbital parameters which give 
aatellite coordinates in an inertial system. The 
determination of these parameters is quite complex and is 
omitted from this thesis. Readers desiring this information 
are referred to Reference 7, 

The defined inertial coordinate system is illustrated 
in Figure 5.5. The origin of this system is tne same as 
that of the earth rectangulur coordinate system defined in 
mection 5.8. The Lr axis coincides with the pole of the 
orbital plane. ‘The KX, axis passes through the perigee 
and the Y, axis completes the right hand system. 

Sell SATELLITE COORDINATE DETENTION <== erie 
navigational satellite systems require accurate orbital 
information so that the satellite coordinates at any point 
in time may be calculated. Since it is not possible to 
predict the position of a satellite with respect to time 
to the required accuracy for more than a day it is impossible 


9 


[3 
to include a satellite ephemeris in a nautical almanac. 
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FIGUR" 5.9 - INERTIAL COORDINATY SYSTEM 
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To overcome this drawback the satellite contains an ephemeris 
weome SvOred in a Mepnetic memory Which is periodically 
updated with new voarameters from the ground. 

The memory is divided into two parts. The first or 
fixed part contains the varameters required to define the 
best ellivtical apvroximation to the actual orbit. These 
narameters cnange slowly with time and sre taken as constant 
over a twelve to sixteen hour period. The descrenancies 
between this ordit and the actual orbit are smail. The 
second or ephemeral part of the memory contains the small 
werrecctions to the hasic orbit for veach time voint. ~The 
orbital parameters stored by the satellite memory are given 

fat 
in Table 3.1. From this data the navigator is adle to compute 
®ie position of the satellite in @mertial wmectanguier 
coordinates at the beginning and end of each time interval. 


ar S 


, ° ° Q 
Yo calculste tne satellite coordinates Avs Via,” (OC ae 
a 


1 


the navigator goes sequentially tnrough the following 
[3] 


calculations where t, is the mean satellite vass time. — 
(ou LG 

Wp = Op 4) p OP 
O t 

JINN = Ax! A yA m 


= VOD (n- Ant We Ah 
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PACBLY “oel ~ See LIE CPRITAL PARAV ST ERS 


Fixed parameters: 
t x time of first perigee after 02 or rehire 
* mean angular motion « 277/period. 
WO ; 

(0) 
O *« argument of perigee at to 
LJ * precession rate of perigee. 
2. 8 CCenericCluy Gf sero u. 

A % mean semi-major avis, 

(o) ° a s ~ 3 
/ \ + ripnt ascension of ascending node at ty. 
‘dy * precession rate of ascending node. 
C; % Cos¢ »%s intlinavion of crbival oie. 
Sec! Oe 


A\M % change of mean anomaly for 1 hour. 


SM « Change of mean anomaly for 2 minutes, 
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™phemeral parameters: 


cr 
f-- 


al. 
bd eal 


time after intepral half hour 


=] 
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3% 


correction to true anomaly at t;. 


Oo, Ory 
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f- 


correction to semi-major axis at Ce 
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Then the satellite coordinates in the inertial system at 


time t; are piven by: 


x5 = AA COs iy a -_) 


i 
+ aa 7 : “7 2 3 | 

(3.44) ve = a, [sin 3,(1 - ©2)2 
Z= —" 6 


These are the satellite coordinates in the defined inertisl 
Seoordinate system at the jo wh time point. sine ercer tome 
macasloOr DOSibion determination in (5S 7fe je tteeem oreo 
ordinates must be transformed to the defined earth 
rectangular coordinate system. 

Oele TRANSFORMATION OP SUNT PAR COC DN 0S eee 
transformation of the inertisl sutellite coordinates to the 
earth rectangular coordinate system can be considered in 
two distinct steos. The first steov involves two rotations 


AnlGees— axea invo vUhne Gaquscorlae 


to Dring the inertial Xt I 


o5 





r Eas COMCOinNcClLae Wom coe Prova loOns lwo 1S 


weet Pook hcscumeitatlons.define the cs 


axes illustrated in Figure 3.5. The second step involves 


plane and the Z 


vee eugiel 2A) 

i ae 
Peiacing the ie axis to coincide with the X axis. This is 
weomelicaimed by the earth's rovution edvout its axis. Since 
the earth rectangular coordinate system rotutes with respect 
to the inertial coordinate system the second step t#akes the 
form of a rotation of the inertial prime system about the 

2 axis through an angle of a = GAST.. Upom comm bepred 
Sievunese steps the satellite coordinates will be in the 
earth rectangular coordinate system used in the position 


Gerivation. The matrix form of the required rotations is 


given below: 


Step 1 - rotation about the 4; axis through an angle 
(0) 
of - Oe followed by. a rotation about the 


! ° 
A+ axis through an angle of -{ . 


(5.51) 


4 


~ 


°) (0) ss Tl eaten 


| ( 
x cos(- U,, ) sin(- Y, ) OF | Xe 


ane 
S 
O 
_ 


- 
: (0) (0) 


© 
O 
TA 
-_— 
i 
oe 
~~ 
na eS 


| 

| 
YF oO er Gey |-sin(-Q- cos(- 7 ) of Yy 

| | 
2 \ -sin(-() je O O 1 ka 


On 
e>) 





Step 2 - rotution about the Z axis through an angle 


of (- S\_- = GAST} ° 





oe ) 
x | cos(/\., - GAST) -sin( A, - GAST) QO | | x! 
N N wie. 
bo: a 
val ee sin( A. - GAST) cos( A» - GAST) Oa ha 
a 
_ 
1 
wf O O i = 7, 
Le ees 
3.15 SUMVNARY - - The theoretical develonment of 


eanations for position determination given in the vrecesding 
focecions illustrates how™aA navisavor ws able to fix hig 
position on the surface of the ellinsoid used as an assumed 
figure of the earth from observations of a navigational 


satellite. Initial assumed coordinates X,, V_; 


anici 

12 
are required from some independent source. In addition, 
the navipator must determine the integrated cycle count 
trom at least three satellite ooOsscrvations ang bees ole to 
SCarculate the C@erdinateawor the sovelixve sue eo pmo o. 
these observations. The ™soluvion volves tne ees orn nae Ter 


Of mccmree lone aX» am and qZ, which sre applied to the 


nN’? 
assumed initial coordinates. This gives the new position 
coordinates in a rectsnpular system. These coordinates 


are then expressed in terme of peodetic latitude and 


Llonmertude. 
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4. OPERATIONAL SYSTEM COMPCNENTS 


4.1 GENERAL. - ~ The Navy Navigation Satellite 
System has been operational since July 1964. It consists 
of three satellites, four tracking stations, an operations/ 
computer center, two injection facilities, and numerous 
navigational receiving units. The ground network is based 
entirely within the United States and is overated oy the 
U. S. Navy Astronautics Group. The locations of the 
components of the ground network are shown in Fifure ait" 

The task of the sastronautics Group is to determine, 
meocesa, and load orbital and time information into the 
memories of the operational satellites. These operations 
must be carried out within critical time limits. Hach 
satellite is tracked every time it comes within range of 
one of the four tracking stations. The data received 
from the satellite is then sent to the computer which 
updates the orbit, predicts the satellite positions for 
the next 16 hours, and generates the orbital parameters 
in a coded message. This message must be sent to the 
injection cAeeieey in sufficient time to be injected into 
the satellite before it has used up all of the information 
previously stored in its memory. ach cycle requires 
about 12 hours per satellite. Since the satellite memory 
igs Limited to 16 hours of information, operations must 


meet an exacting schedule. 
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4.2 ‘VHH SATSLLITE. - - As was previously mentioned, 
details concerning the operational satellite are not readily 
available due to security restrictions. However, it can be 
assumed that the oversxtional satellite employs the basic 
characteristics develonved by the Transit research and 
develooment series described in Ghapter @. 

4.3 TH TRACKING STATION. - = The four overational 
tracking sttions sare located in Maine, Minnesota, California, 
ema Hawaii. The tracking station seunuipment consiats of a 
Westinghouse AN/BRN-3 type receiver, a doppler cycle counter, 
an ultru-stable oscillator, «and communications equipment. 

The data obtained by the tracking station is put into standard 
teletype format consisting of two parts. The first part is 
called the header and contains the identifying information 
reouired vy the computer. The second part contéins the 
actual data. The data consists of a series of ll-digit 
data points. ach point represents one doppler reading and 
the time that the reading was taken. “hese readings are 
usually taken every four seconds. Generally 100 or more 
data points are obtained during the pass of one setellite 
within. radio. line-of-sipnt sot the Gieeei 1 pests Bien aD 
messages gent to the operations/computer center, the first 
five digits of a dita point represent the time of day in 
seconds (UT) and the last six digits represent the number 


of microseconds it takes to receive a predetermined number 
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of cycles from the satellite. This preselected value is 
such that the time aporoaches but never equals one second. 

In order to send the great quantity of data gathered 
by the tracking stations to the overations/computer center, 
high speed communication equinment is required. A repid 
data-link called TRADAT, which has a transmission rate 
five times faster than teletype, is presently being used. 
TRADAT is essentially error free because it incorporates 
self-checking circuitry. 

4.4 OPZRATIONS/COMPUTSR CENTGR. - - The operations/ 
computer center is located at the headquarters of the 
Ue S. Navy Astronautics Group at Point Mugu, California. 
All deta from the tracking stations and the computer pass 
through the operations center. This center at all times 
maintains the operational system status and the positions 
of the satellites as well ag related sy:tems information. 
It is the overations center's responsibility to insure 
that the tracking stations track each satellite and that 
sufficient data is provided to the computer to allow 
accurate orbit determirations and predictions. It schedules 
the computer center runs and insures that the orbital and 
time parameters are gent to the injection facility in time 
to meet the injection pass of a satellite. 

Four satellites in orbit will produce an average of 
84 messages ver day from the tracking stutions. The data 


in these messages pass through the operaticns center and 
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are sent to the computer center whose Vriidenmy so0 1s 
S2oe UDdeting “nd position prediction.* Pme computer 
center has an IBM 7094 II/1410 confipureation with a 
360/40 for commend and control functions. The dovpler 
data are fed into the computer which is programmed for 
orbital analysis. The machine computes a theoretical 
orbit of the satellite and compares a theoretical Gcopler 
Curve obtained from this orbit with the actual domsiler 
curve. The residuals are squared and summed and 2 new 
set of corrected variables is obtained. These aliow the 
computation of a new theoretical orbit anc & new tnhsorevica, 
doppler curve to be again compared with the observed dati, 
Tnois process is repeated until the sum of the scausres of 
Gao csTauals 1s a Mintrur 20 chis DoLrt™isc 1a escunmec 
that the best determination of the oroit has been made. 
In this way orbits of high precision can ode determined 
and the future voeition of the satellite can be or3idicte? 
Wilt" the required accuracyd?! It takes “nearly two hours 

j 
GO ucdate a sineplie satellite ornit and prepare an injection 
message. With each satellite requiring this informution 
aporoximately every twelve hours, the computer wili be used 
to capacity. ‘The messages prenvared by tne computer contains 
24,960" Dits of crbhivaieeand Cine =n om bionm Waren comer ce 
the message proper and 220,000 bits of information which 
control the eavuismenv av the Giayecouron oct) 1 Uy semen 
satellite pass. 
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‘Although TRADAT speed and ccuracy are ourticulnarly 
adé9oted to communications odetween the tracking stations 
and the operations/computer center, it still does rot 


rovide the speed and accuracy levei reauired to handle 


3 


he injection messages which are sent from tne opnerations/ 


ct 


computer certer to the injection facility. These must be 
HOO per cent error tred. High-speed eavipment with a data 
rave more than twice that of TRADAT is used fer the sole 
purpose of sending injection messages, which, from the 

time tney are fed into the equipment, are entirsly sutomated 
until they are transmitted by tne injection facility 
antenna to tne satellite. 

4.0 THe INJSCTION FACILITY. - - The operational 
system employs two injection facilities, one in Californic 
soemoenc in Minnesota. Hacheinjeetion [ee1 1s 2yecemoawa. 
eee 1000 Sesmcolic directional antenna in conjunc lion 
with a 10 kilowatt transmitter. A small computer similar 
to that used in the AN/B°N-3 submarine navigaticn éauinment 
Pentrels one operation. Thew™imjection Wer emac Tre asthe 


operutions/comouter center is recelvac ev the injecvion 


} 


Facility and stored. ‘“yhen the sctel.ite is in range of tne 
mijection facility the messape proper Containing O.td or 


Earormation in combitm@eanns of Dbinaey em] ame ros is 


ti 


re 
Converten Co WavVetorm co modulate the Gurricr WAVe Of Une 


transmitter. Tne message is transmitted into the satellics 
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in 15 seconds. Since the satellite and the earth are in 
motion during transmission, the distance from ths transmitter 
to the satellite varies constantly. Therefore, the time 
reouired for the signal to reach the satellite must de 
comcuted in advance and the bit rats must be adjusted to 
ensure that tne satellite receives them at s 2ven rate. 

The bit rate is changed 32 times during the 15 second 
transmission. The injection facility also reads the message 
as retransmitted from the sutsllite to determine the success 
of the injection attempt. If necessary, the injection 
facility recycles to do the wnole opneracicn over again. 
complete cyCie Tor eacn injection requires two minutes. 


- 


The egquinom-nts of the injection facility consist o 


4) 


commercially available components wherever possible. = 
fact the transmitter is a modified commercial televisicn 
Baittor 

4.6 NAVIGATION “SQOUITPMANT. - - Two types oF shipboard 
navigation equipments Nave been developed, ones for submarines 
and ths other ror susctace ships. Tae AN/B-S was deveidnedc 
by Westinghouse and is designed specifically for suvubmarin: 
use. It incorsoratves many features adapted to submarine 
environment, and certain pieces of inveri:ce equipment to 
Perms UsoM eel Lonmin CONnIMAClL OMe Bem ove 1 OT: 
equioments presently used abourd Fleet Bu:listic pissile 
submarines. The AN/BRN-3 is completely sutomatic in 
Operation. It comoutes it's own slerts of wren navigational 
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satellites will be witnin radio range of the suvomarine 
position. It then computes the satellite position for 
specific intervals during a pass based upon knowledge it 
already has of the sateliite crbit. Using an assumed 
submarine position it precomputes the expected donpier 
curve. When the actual doopler curve is received the 
AN/BRK-3 adjusts and readjusts the assumed position until 
the closest possible matcning of the two curves is achieved. 
Since the assumed position is supplied by an Inertisl or 
some other high quality positioning system, the curve- 
Fitting process takes olsace almost in real time. The 
position at which tne two curves fit most closely togethor 
momen to oe the submarine's position. “Whis fix as 


dle. 


transmitted to tne other navigation equipinent on board, 


a) 
{2 


necessary. In addition, the esoulpment performs its own 
readiness checks, diagnoses any trouole, und prints out in 
english what part to replace. 

The equipment developed for surface ships ie designated 
the AN/SRN-9. ic is far less complex than the suomarine 
aquipment since it does not have the environmental znd 
interface oronlems associsted with undersea operntions. 
The AN/SRN-9 receives the signs: from tne satellite and 
sutomabtically computes tne ship nosition. It uses the 
orbital parameters provided by the satellite to compute 
the part or the orbit correspond inigvo wre eon! Coeur 
received. From this it computes tne satellite coordinates 
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aCethe observation time points. with this information, 
olus the integrated doppler cycle count for at leust three 
2 minute time intervals, it computes the corrections dx» 
ay: and dA, to the assumed position coordinates. The new 


position coordinates are then used to calcuiate the navigator 


Peaeetuae und longitude. 
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>.» SYSTSN ACCURACY 


5.1 GSANSRAL. - - In the course of the development 
of the Navy Navigation Satellite system, individunl error 
components were anslyzed in order to determine and improve 
those components which most affect the system. Thtoretical 
values based on assumed ideslized conditions were determined 
for most components. These values can be expected to vary 
under overating conditions. Specific values pertaining 
to tne accuracy of fix determinations based on the use of 
the prototype and operational systems are not aveailsble 
Gue to security restrictions. However, several tleet 
units have made extensive use of the operationai system 
and fix accuracies of t 0.1 mile or better have been made 
Pool iC. 

From the navigator's point of view the primary error 
components of the operational system ares: 

1. Yonosnheric refraction error. 
oe Freauency error. 
Ope Limemerwors 


4. Grror in the determination of the velocity 
of electromagnetic propagation. 


Oo. Navigstor instrument error. 
6. Navigator velocity error. 
7. Satellite coordinate error. 


A short discussion of these error components is given in 


Oe! 





the following sections. It snould be noted tnst due to 

the nature of the system, the satellite coordinate error 
Gan be oroken down into individual error components similar 
to these encountered in fix determination. The principsl 
difference is in the high order gravity coefficients 
required for orbit determination and prediction. 

9-2 ILONOSPHGRIC RuFRACTION. - = In Section 5.2 it 
was shown that a first order correction for ionospheric 
refraction can be determined by the transmission of two 
coherent frequencies from the satellite. The resulting 
doppler shift still contains higher order effects of the 
series (1/f.)- In order to determine whether contributions 
higher than the first order are significant, the theoretical 
effect of refraction on the doppler shift was consincerseé 
in Reference 9. Table 5.1 summarizes the results pertaining 
to upper bounds on the various contributions to the refracted 
doppler shift for various frequencies within the transmissicn 
ee xcept for the vacuum term, typical values are 
premsoly a factor of three wo Cerniwemalicre ‘ie Lirguaorder 
term is proportional to the time rate of change of the 
total number of electrons in 4a unit cross section tube 
extending from the observution point to tne satellite. This 
is the term that is eliminated when two frequency doppler 
Gata is used. The second order term considers Faraday 


rotation and its sign denends on whether the antenna is left 
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Or Might circulary polurizod. The third order correction 
is composed of four terms. The first term arises from 

the fourth order term in the expvansion of the refractive 
index series. The second term revresents the difference 
Detween the optical path and the geometric path. The last 
two terms reoresent modifications to the solution of the 
wave equation using the geometricsl optics approximation. 
These terms depend upon vuricus derivatives of the electron 
density. Unless the satellite frequencies are below VHF 
the two terms not obtainable from the geometrical optics 
approximation are negligible. 

Table 5.2 presents the estimated maxims of tne second 
and third order contributions to the vacuum doppler shift 
when dual frecuency dita are combined to eliminate the 

[S| 
first order refraction term. Three typical frequency 
combinations srs given and the table entries have heen 
converted from frecuency shift to errors in the runge rute 
in meters per second. 

Examination oF the tables indicates that if the lowar 
of the two frequencies is not less than 100 mc/s and the 


Llonosphere is not extremely disturbed, tne second and third 


order’ rotracthonscontwabubdons eaewkdebeemecrlirible. ° 
Oed FRESUsNCY AND Tie “GRRORS. = = Whe navigavor 
must have an accursute clock and 4 reference freauency in 


order to determine nis position from doppler observations. 
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In addition he must know the frequency transmitted by the 
Sewollite and the satellits position 4s a functicn Ge Gime. 
In Chapter 3 it was shown that to the required accur.cy the 
frequency of the reference oscillator is assume? unknown 
and the position determination computations were arranged 
20 as to eliminate tne need for this value. The navigator 
is able to recover time from each satellite observation to 
an accuracy of better than one millisecond. This time is 
maintained by & clock derived from the reference oscillator 
whicn has a drift rate of less than 1 pert in to per day. 
Therefore, the time and freouency errors in the Navy 


Navigation Satellite System both derive from the drift of the 


satellite oscillator. <A drifting oacillator g 


ae 


Vea ores 


s in error 


|e 


track navigation error beceuse the frequency 


J 


em €o2uses a time error in the clock that it is driving. 


This time error ?lso gives an slong trick error. Navignticn 
epmors Of about 1 kilometer"result from an error of | oar 

eye 10° in the satellite frequency. However, satellite 
freaquency can be determined with an iccurscy of at Least 


10 


Meearcein 10 and predicted with un accurtcy of at least 


* O hd e J 
feos in ioe per day. Since the sstellite oscills«tor 
can be tuned by command from an injection staticn twice 
per day, it's frenuency can be kant to a predetermined 

° . 4 bee s @ m™m 
value within an accuracy of a least @ perts in 10°. This 


[io] 


causes a navigation error of poroximsately 20 meters. 
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54 VeLOCITY OF 2LaOTROMAGNETIC PROPAGATION. - - The 
value for the velocity of slectromagnetic enersy vropaugation 
memoalcon FCO be a constunt in enlenlatio@me Por nosi tion 
determinetion using the Mavy Wavigation Satellite Syetem. 
The value is probably known to an accuracy of 1 part in 
e000 ,000. This would vermit position determimation to an 


rf 
’ 


accuracy of «about two meters. 


a4 


Future improvements in this 
constant can be accomodnted by modification of the system 
based on the current accepted value. 

Oe) NAVIGATION INSTRUMENT GRROR. - - In order to test 
the contriodution that the navignution equioment rakes to the 
navigation error, observations were made on experinental 
satellites transmitting on four frequencies. Two névigation 
instruments at the same site mace simultaneous observations 
on the satellite. One ineatrument used two of the satellite 


Freaquencies. 


frequencies and the other used the remuining two 
The difference hetween the navipsited vositions is then tre 
result of two inoependent measurements of vosition. The 


a 


ite 


Mmesutt.of one setof tests obtained with tie Tea 


satellite gave a mean difference in position of 50 meters 


(3 ; 

jm one coordinate and, 360 meters in the other inc ag 
axpected that inetrumentution hse imnroved since this time. 
5.6 NAVIGATOR V=LOCITY GRRORs« - = Simvee the Wawi pager 


must use u measured value of his velocity to determine his 


position at the times of the satellite obdgervations, an 
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error in this measmirement will cause sn error in the fix. A 
ship acted on by unknown current, windapge drift, send ses 
condition effects amounting to a one knot velocity error in 
the northerly direction has an estimated position error of 
12 

O.5 nautical miles. A one knot easterly velocity error 
causes a nepligible positional error st mid-latitudes since 
the satellite is in a polar orbit. It was determined to at 
least three place accuracy that the position errors varied 
linearly with the velocity ue This does not appesr 
to be a gevere limitation for surface snips and submarines 
where accurate pit-log or inertial system velocities are 
normally available. 

oof SATSLLIT@ COORNDINATS SRROR. - - Apart from the 
preceeding, error comnonents which are normally nepligibie, 
the only remsining source of system error is in the satellite 
coordinates transmitted to the user. Any error in thesé 
coordinates is reflected almost exactly one for one in the 
eseeuracy of Che navigation fix. Tests onthe eize of sense 
errors in the predicted orbits of the Navy Navigation 
Satellite System were conducted by computing and predicting 
the orbits at four tracking stations. These predictions 
were made using geodetic parumeters in use at the Naval 
Weapons Laboratory during four periods of time tnrough 1964. 
The predicted orbit wes comptired with an orbit based upon 


independent doopler ohservations mide by world-wide trackin,y 


stations during the prediction interval and using the latest 
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aveilable pravity coefficients. The resulting tweive 
hour prediction errors for random sets of observations 
made on polfr orbiting sextellites show a decrease in 

size from over 1000 meters in 1960 to the value of 50 


A 
mevers in 1965. 
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9e8 PROPCSHD “RROR ANALYSIS. - - A high speed comouter 
may De easily programed to analyze the effect of indivicual 
and/or collective error components on a given satellite 
position determination. The hasic inouts of this program 
would consist of the following porameterss: 
Le ree ee 


n var 
2. 4iax, LdAy, itz 


n 


me AR 
‘ia iP 

SRM? 2 

Ge. C 

eee 


The satellite coordinates as a function of time in an 

earth centered rectsngul‘ir coordinate system are availabdle 
from a Geodetic Science Computer Librery Program at the 

Ohio State University or from some other similar program. 
Wsing satellite coordinutes bracketing three or more time 
intervals and choosing the above parameters which ere assumed 
errorless, we can determine the expected integrated cycle 
count from (3.7). Assuming errorless instrumental and 
environmental conditions the calculated cycle counts and 

the above parnmeters may be considered as an ideal errorless 
vosition determination based on satellite observations. 

This can be verified by substitution of these values into 


the eauations developed in Chapter 5. The position 
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coordinates voroduced by calculation will be equal to tne 
initial navigntor coordinates since the corrections ax 
qy,» and dAf will be zero in our assumed errorless 
condition. We have thus created an ideal errorless 

position determination which can be used as a standard 

in an error anvlysis. By random and/or systematic 
alteration of the assumed errorless input parameters we 

can determine their individual and/or combined effect on 
the final position determination. A study could also be 
made to determine the effect of decreasing the time interval 


and using more observations in a least squares solution. 
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Ge. SUMMARY 


The adaption of spuce technolory vrovides a new apvroach 
to the earth-bound problem of marine navigation. The Wavy 
Navigetion Satellite System is based on the donvvler shift 
in 6 radio frequency transmitted from & near-earth satellite. 
By exact measurement of the doppler shift it is possible to 
determine the location of the receiver on the surface of tne 
assumed figure of the earth. The accuracy obtained by using 
this technique is possible because the satellite is very 
precisely restricted to an orbit defined by physical laws 
of motion and the quantities observed, time and frecuency, 


10. Of -oaal 


can be measured to an uccuracy of one part in 10 
the possible satellite paths permitted there is only one 
which results in a particular curve of doppler shift. 

Two basic vroblems are inherent in this system. One 
is that the ionosphere bends radio waves and thus fives a 
false position of the satellite. This hss been overcome 
by the use of two satellite transmission freouencies whicn 
allow a first order correction to be comnuted. The second 
basic problem is the determination of the pravity field of 
the earth. This affects the determination and prediction 
of satellite orbits. Tnis problem has been recuced to an 
acceotable level by using satellite orbital data to calculite 
new gravity coefficients. 


The use of a navigitional satellite system based on 
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the dopnler shift of VHF radio weves appears to nrovide 

the best present day solution to the military navigation 
recuirements. World-wide coverage ie obtained vy use of 
ool1r orbits. The system is inherently all weather and 
immune to interference since line-of-sipht radio frequencies 
are used for both tracking and navigation. Since the 
satellites transmit their ephemerides the navigator need 
neither interrogate the satellite nor receive orbital and 
time information by other communication links. ‘The system 
igs also reliable because a temporary inverruption of 


reception during 4 pass does not preclude the use of 


fix. It is assumed that the system meets the milivcary 


accuracy reauirements. 
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